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We briefly review the calculational procedure for the PQCD prediction for hard exclusive quan- 
tities and reconsider the problem of the factorization scale dependence. 
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I. INTRODUCTION 

The application of perturbative QCD (PQCD) to ex- 
clusive processes at large momentum transfer is based on 
the factorization theorems P, L3- EI ■ The main idea is 
the separation of short- from long-distance effects in a 
sense that the high-energy region, being highly off-shell, 
is factorized from the low-energy region, which is charac- 
teristic of the bound-state formation. The factorization 
may be carried out order by order in perturbation theory. 
The information concerning the long-distance dynamics 
is accumulated in the distribution amplitude (DA), one 
for each hadron involved, whereas the short-distance dy- 
namics is represented by the hard-scattering amplitude. 
The separation of the short- from the long-distance part 
occurs at the factorization scale which is usually chosen 
by convenience. Furthermore, PQCD calculation to the 
finite order necessarily requires the renormalization of 
the UV divergences and introduces therefore the renor- 
malization scale dependence in the final result. 

As is well known, one of the most critical problems 
in making reliable PQCD predictions for exclusive pro- 
cesses at large momentum transfer is how to deal with 
the dependence of the corresponding truncated pertur- 
bation series on the choice of the scheme for the QCD 
running coupling constant as(/z fl ) and on the choice of 
the renormalization scale /zr, as well as, the factoriza- 
tion scale (If. Although the physical quantities depend 
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neither on the renornralization nor on the factorization 
scale, the PQCD prediction at the finite order bears the 
residual dependence on the renormalization and factor- 
ization scales, the choice of which introduces theoretical 
uncertainties in the prediction. 

A lot of work has been devoted to the analysis of the 
renormalization scale and scheme dependence [a, la, 13 • 
The problem of finding the optimal renormalization scale 
in a given scheme has been widely discussed in the liter- 
ature and, apart from the pragmatical choice n R equal 
the characteristic scale of the process, three quite dif- 
ferent approaches have been proposed: the principle of 
fastest apparent convergence (FAC) Q, the principle of 
minimal sensitivity (PMS) [g and the Brodsky-Lepage- 
Mackenzie (BLM) scale setting \j\. A physically moti- 
vated formalism in which any two perturbatively calcu- 
lable observables can be related to each other without 
any renormalization scale or scheme ambiguity has been 
developed || 

In contrast to the renormalization scale, somewhat less 
attention has been paid to the role played by the factor- 
ization scale. Although one can encounter in the liter- 
ature several extensions of the treatments of the renor- 
malization scale to the treatments of the factorization 
scale H, when examining the hard exclusive processes a 
convenient choice /ip equal the characteristic scale of the 
process, i.e., the large momentum transfer denoted by, 
say Q 2 is mainly used, with the justification that for such 
a choice ln(Q 2 / /j, F ) logarithms, giving rise to the growth 
of the coefficients in the expansion of the hard-scattering 
amplitude when Q 2 ^> n F , vanish. Obviously, the result 
will differ for some other choice of fip. Then one can try, 
similarly to the renormalization scale problem, to justify 



other choices for the factorization scale by examining the 
underlying dynamics of the process |l(J, LLLi ■ 

In this paper we review the prescription and the in- 
gredients of the higher-order calculation of the hard ex- 
clusive quantities (obtained in the, so called, standard 
hard-scattering approach |TJ.yy. 3) and we reexamine 
their factorization scale dependence. We show that the 
residual factorization scale dependence in the finite order 
of PQCD calculation reflects the failure of the proper re- 
summation of all hi(iip) logarithms. Thus, taking into 
account the factorization scale dependence of the hard- 
scattering amplitude and of the distribution amplitude 
by consistently including all terms that are effectively of 
the same order in as, the PQCD prediction for an exclu- 
sive one-scale process is free of any residual dependence 
on the factorization scale at every order of the PQCD 
calculation. The unavoidable theoretical uncertainty of 
a particular order of the PQCD calculation remains to be 
only due to the renormalization procedure. Nevertheless, 
we comment the problems one is left with when adopt- 
ing such a procedure, especially in the case of multi-scale 
processes. 

In Sec. [H] we introduce the ingredients of the stan- 
dard hard-scattering picture on the example of the pion 
transition form factor, while in Sec. Illll the higher-order 
calculational procedure is outlined. The discussion of the 
factorization scale dependence is given in Sec. IIVI Sec- 
tionals devoted to concluding remarks. 



where ijl 2 r is the renormalization scale. The dependence 
of the coefficients of the expansion J5J on the scales n R 
and p, F is of the form ln™(^/<2 2 ) and In™ (p 2 F / Q 2 ) , re- 
spectively. 

The pion distribution amplitude $(#, lip), although in- 
trinsically nonperturbative, satisfies the Brodsky-Lepage 
(BL) evolution equation 

d 

(4"frT®( x > &) = v ( x > u > Ml) ® $(«, a4) » ( 3 ) 
o[i F 

where V(x, u, fip) is the perturbatively calculable evolu- 
tion kernel 



V(x,u,/J, F ) 
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V 2 (x,u) 



(4) 
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The solution of Eq. J3J can be represented as 

(5) 
where & LO and & NLO denote the leading order (LO) and 
next-to-leading order (NLO) parts, respectively. When 
convoluting the finite-order results J2J and (J3J) according 
to (JTJ, one is usually left with the residual dependence 
on both fift and lip. The origin of the latter will be 
explained in the following. 



II. STANDARD HARD-SCATTERING 
PICTURE AT HIGHER-ORDERS 



For definiteness, notational simplicity, and clarity of 
presentation, we consider the high-energy behavior of the 
simplest exclusive quantity, the pion transition form fac- 
tor F 7 k(Q 2 ), defined in terms of the J*(q, fi) + j(k,u) — > 
7r(P) amplitude. For large momentum transfer Q 2 {= 
l ), the general factorization formula (U- LJ ll Q| for 



F-yniQ ) reads 



F^{Q 2 ) = $*(z,/4) ® T H (x,Q 2 ,fi 2 F ) 



(1) 



Here, $>(x,iip) is the pion distribution amplitude; 
Th{x, Q 2 ,iip) is the hard-scattering amplitude; [ip is the 
factorization scale, and x denotes the pion constituent's 
momentum fraction, while ® = L dx. 

The hard-scattering amplitude (HSA) Th can be ex- 
plicitly calculated in perturbation theory and represented 
as a series in the QCD running coupling constant as(fJ, R ) 
by 



(0)/ 



T H {x,Q\v, F ) = T%>{x,Q i ) + 



Q g(/4) rp(l) 



4tt 



Tg'(x,Q',n F ) 



»I(m|) ^(2) 

(4ir) 2 



± H (x,Q iUf'^r) 



(2) 



III. CALCULATIONAL PROCEDURE 

In order to be able to examine the origin of the resid- 
ual dependence on fip, we first reexamine the calcu- 
lational procedure and the ingredients of the standard 
hard-scattering picture for F 17T (Q 2 ). 

The HSA Th is obtained by evaluating the 7* +7 — > qq 
amplitude, which we denote by T. Owing to the fact that 
final-state quarks are taken to be massless and on-shell, 
the amplitude contains collinear singularities. Since Th 
is a finite quantity by definition, collinear singularities 
have to be subtracted. Therefore, T factorizes as 



T(u,Q 2 ) =T H (x,Q 2 ,fi 2 F ) ® Z T ,coi{x,u\^ 2 F ) 



(6) 



with collinear singularities being subtracted at the scale 
H F and absorbed into the constant Zr.coi- The UV 
singularities are removed by the renormalization of the 
fields and by the coupling-constant renormalization at 
the (renormalization) scale fi R . 

The process-independent pion DA in a frame where 
P+ = P° + P 3 = 1, P~ = P° - P 3 = 0, and P x = is 
defined [HHEI as 



$(u) 



dz 

~27T 

x (0 



i(u-(l-u))z-/2 



*(- z )Zj*n*(z) 



■ (7) 
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where $7 = exp < ig J_ ± dsA + (zs)z /2 > is a path-ordered 

factor making $ gauge invariant. Owing to the light-cone 
singularity at z 2 = [j, [jj| the matrix element in Q is 
UV divergent. After regularization and renormalization 
at the scale jl 2 R: z 2 is effectively smeared over a region 
of order z 2 = — z\ ~ l//&fj- As a result, a finite quan- 
tity, namely, the pion DA $(i;,/2|j), is obtained and cor- 
responds to the pion wave function integrated over the 
pion intrinsic transverse momentum up to the scale jl 2 R . 
The pion DA as given in Q, with \w) being the physical 
pion state, cannot be determined by perturbation theory. 
If the meson state \ir) is replaced by a \qq;t) state com- 
posed of a free (collinear, massless, and on-shell) quark 
and antiquark (carrying momenta tP and (1 — t)P and 
pseudoscalar meson quantum numbers), then the ampli- 
tude 10 becomes 



,„./, = / ^_ e i(«-(l-«))^"/2 



x (0 



*(-z)X_^fi*(*) 



qq;t 



(8) 



pion bound-state dynamics) determined at the scale /1q, 
while (f>v(v,s;ji R ,(j,Q) governs the evolution of ^(v,/j^) 
to the scale fi 2 R . By differentiating i|ll[l with respect 
to jl R , one obtains ©, with the evolution potential V 
given by 

V{tf R ) = -Z^ ren (ti) f Afl^-^,r«n(AK)) ■ (14) 

To simplify the expressions, the convolution ((g)) is here, 
and where appropriate replaced by the matrix multi- 
plication in x-y space (the unit matrix is defined as 
1 = S(x — y)), while the x, y variables are suppressed. 

By convoluting the amplitudes T(u, Q 2 ) and $(u), 
Eqs. iJBJ and JTTJ, respectively, in analogy with 0,13 
we obtain the pion transition form factor F^ 7T (Q 2 ): 



F JW {Q 2 )=&{u) T(u,Q 2 ) 



(15) 



Taking JSJ into account, we can express Eq. JJJ as 



Now, in order that the factorization holds, p, R has to to 
coincide with fi F and by making use of the fact that 

Zt,coi(x, w, h%) ® Z^ ren (u, v; n 2 F ) = S(x - v) , (16) 



$(w) = 4>(u,t)(g> {qq;t\ic) 



(9) 



The distribution 4>(u, t) can be treated perturbatively, 
which enables us to investigate the high-energy tail of 
the pion DA and its evolution. The (f>(u, t) distribu- 
tion is multiplicatively renormalizable owing to the mul- 
tiplicative renormalizability of the composite operator 
\l/(— z) 7 + 75 ft &(z). This means that the UV singular- 
ities that are not removed by the renormalization of 
the fields and by the coupling-constant renormalization 
factorize in the renormalization constant Z, 



(p,ren 



at the 



(renormalization) scale p, R . Apart from UV singularities, 
the matrix clement in JSJl contains also collinear singular- 
ities. Subtracting these singularities at the scale fi 2 and 
absorbing in Z^ :CO ;, we can write Eq. (JSJ as 

4>(u,t) 
= Z^. ren (u, v; jl 2 R ) ® (j>v{v, s; Jx R , /tg) ® Z$ tCol (s, t; $) . 

(10) 

By combining 10 and IjlUfl . we obtain the distribution 
<£(m) in the form 



$(li) = Z < j, ! re n (u,V,iJ,%)®§(v,p? R ), (11) 



where 



Here. 



r-,2 .,2 



*(»,A*ji) = M«.*;Afl,/*o)®*(*.A*o)- ( 12 ) 



$(s, /ft) = Z , co ,(s, t; /^) $ ( TC ; *|tt> (13) 

represents the nonperturbative input (containing 
collinear singularities and all effects of confinement and 



the divergences of T(u, Q 2 ) and $(u) in (|15|) cancel (this 
has been explicitly shown in | 15| up to n/ proportional 
terms of 0(a 2 )) and we are left with the finite perturba- 
tive expression for the pion transition form factor Jy. 

We note here that the same factorization (and renor- 
malization) scheme is employed in the hard-scattering 
and DA part, i.e., in Eqs. |JB} and i|TT|) . respectively. 
Furthermore, as pointed out in [l6j, the evolution equa- 
tion as defined by J3J and (0} corresponds to the simpli- 
fied scheme fixed by the preference that the distribution 
amplitude should have no dependence on the renormal- 
ization scale 1 . 

It is worth pointing out that the scale n F representing 
the boundary between the low- and high-energy parts in 
||TJ plays the role of the separation scale for collinear sin- 
gularities in T(u, Q 2 ), on one hand, and of the renormal- 
ization scale for UV singularities appearing in the per- 
turbatively calculable part of the distribution amplitude 
<&(u), on the other hand. 

The calculational procedure explained above is illus- 
trated in Fig. n] 



1 Note that, in general, such a residual dependence appears along 
with the evolution kernel depending on two scales: 
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Vi(x,u)+ a °;^' \V2(x,u) - l3 Vi(x,u)ln( ^ 



Here fi% corresponds to the scale of the coupling constant, while 
p? R = fip denotes the scale at which remaining UV divergences, 
due to the renormalization of the composite operator, factorize. 




(l-u)P 




FIG. 1: Pictorial representation of the calculational ingre- 
dients of the pion transition form factor. T represents the 
perturbatively calculable 7* +7 — » qq amplitude , while $ de- 
notes the (unrenormalized) pion distribution amplitude given 
by (|7f, which can be expressed, as in ©, in terms of the 
perturbatively calculable part (j> ^ and the perturbatively 
incalculable part. 



IV. FACTORIZATION SCALE DEPENDENCE 

We next turn to the discussion of the /sp dependence 
of the pion transition form factor defined as in (fT|). 

Concerning the pion distribution amplitude $(#, H F ), 
its dependence on fip is specified by the evolution equa- 
tion © and, as can be seen from Eq. (|12J) . this depen- 
dence is completely contained in the evolutional part <fiv ■ 
By calculating the perturbatively obtainable amplitude 
4> [15| directly from (JSJ, the result obtained for cj)y can 
be organized as 



<fo (Mf>Mo) = ^(mI.Mo) 



osJPf) a nlo 



(Mf.Mo) 

(17) 



where 



iLO 



F(Mf,Mo) = 1 



asdufO 



+ 



47T 

«|(Mf) 



ln4^i 

Mo 



iNLOr 2 2\ 

</V (Mf,Mo) 



(4tt) 2 

+ ••• , 

as i^ 2 



In 2 ^ i (V? 

Mo z 



T^ln4^2 
4tt ^ 2 



•&Vi) 

(18a) 

(18b) 



and functions V n represent the n-loop evolutional kernels 
appearing in Eq. Q. The terms explicitly given in 118() 
correspond to the results of the two- loop calculation [15( . 
In writing Q18[) . the use has been made of 



a s (n 



4tt 



Mo A> 



"s(m 



as(^o) 



= 0(a° s ) 



be determined by solving the evolution equation (JSJ or 
equivalently 

2 J. / 2 2\ 

Mf-5-^-<M u , s ,Mf>Mo) 

= V(v, s', mf) O 0v(*', s, mI, Mo) • ( 19 ) 
The LO result is of the form 



<^>, S ;Mf,Mo) = E'^^^ 72 ^- 1 ) 



xCy 2 (2 s -l) 



((xs{Hf) 



\as(v 2 ) 



(20) 



where N n = (n+l)(ra+2)/(4(2n+3)), while C*^ /2 (2x-l) 
are the Gegenbauer polynomials ( the eigenfunctions of 
the LO kernel V\ with the corresponding eigenvalues 7™ 
|Tl|). The above given complete LO prediction repre- 
sents the summation of all (as In Mf / Mo)™ terms from 
The complete formal solution of the NLO evo- 



lution equation was obtained in (17| by using conformal 
constraints and the form of 4>y LO (corresponding to the 
resummation of (llSbfl ) can be extracted from the results 
listed in [Hi- 
lt is important to realize that the method employed 
above to study the \xp behavior of <fiv can be used to 
examine the dependence of the hard-scattering amplitude 
T H (x,Q 2 ,n 2 F ) on the scale /ip, as well. 

By differentiating QJ with respect to n F and by tak- 
ing into account J2J, one finds that the hard-scattering 
amplitude satisfies the evolution equation 



/4 



8 



T H (x,Q 2 ,[i 2 F ) = -T H (y,Q 2 ,n 2 F ) ® V(y, x; n%) . 

(21) 

This equation 2 is analogous to the DA evolution equa- 
tion J3J . Similarly to the above discussed solution of the 
DA evolution equation, the finite-order solution of H21|) 
would contain the complete dependence on /ip, to given 
order in as, in contrast to the expansion @ truncated at 
the same order and containing unresummed logs. Let us 
note that the explicit expressions for the hard-scattering 
amplitude Tjj (x , Q 2 , fip) in a form J2J), evaluated up to 
^/-proportional NNLO terms, are given in |l5J|. 

The n F dependence of Th{x, Q 2 ,H f ) can be, similarly 
to JT2J, factorized in the function cj>v (y , x , Q 2 , /j, f ) as 

T H (x,Q 2 ,fi 2 p)^T H (y,Q 2 ,fi 2 p^Q 2 )®^ v (y,x,Q 2 ,fi 2 p). 

(22) 
Using (|19fl one can show by partial integration that 122|) 
indeed represents the solution of the evolution equation 

(ED- 



On the other hand, the complete LO and NLO behav- 
ior of 4>v(v, s; n F , Hq) and, consequently, of <&(w, jj, f ) can 



2 The Eq. 1'2H can be also obtained by combining Eq. JSJ with 
Eqs. (HJ and JTC}. 



When calculating to finite order in as, it seems not 
quite consistent to adopt in the literature often encoun- 
tered procedure in which the $>(x,/j,p) distribution ob- 
tained by solving the evolution equation Q is convo- 
luted with Th(x,Q 2 , fip) obtained by the truncation of 
the expansion (J2J. In the latter case, only the partial 
dependence on jip is included (logs are not resummed), 
in contrast to the former case and hence the residual de- 
pendence on the factorization scale [i 2 F enters. 

The proper procedure would be to convolute <& l(T2|) 
and Tjj (|22fl in terms of the same function <j>v, where 
4>v can be given by 118J I with unressummed logs or can 
represent the solution of Ijl9|) . i.e., the resummed result. 
In both cases the n F dependence of <& and Tjj completely 
cancels out and there is no residual dependence on jip. 
One usually uses the resummed form of 4>v in DA $, 
i.e, $ is taken as a solution of evolution equation, and 
therefore this procedure should be applied for Th as well. 

Substituting l(T2|) and i(22|l in Q, we obtain 



= T H (y,Q 2 ,Q 2 ) 
where 



MlM.Q 2 .Mo)® **0>.Mo). (23) 



4>v(y,x,Q 2 ,[i F ) ®0y(x,s,/i|,^o) = 0v(y,s,Q 2 ,/j%) 

(24) 
has been taken into account. It is important to realize 
that the expression H24[> is valid at every order of a PQCD 
calculation 3 , and hence even the finite order prediction 
for F 1 . K {Q 2 ) does not depend on the choice of the /j, f 
scale 4 . 



3 The Eq. 1231 can be easily checked to the NLO ord er IT^I by 
using the LO result |2UJ and the NLO results of Ref. Ilj. 

4 Let us, following Eqs. ||} and 1171 . define the finite order quan- 
tities 
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( 47r )n+l " 



(here n = 0, . . .). The functions <t>ty' (Q 2 , H F ) and $$' (jjl f , $) 
represent the solutions of the evolutional equations 

/4-A^ n) 04."§) = Wife) 9 $»(&,$), 



» F drt vv 



'ftp 
It is now easy to prove that the convolution 

indeed does not depend on fi^,. 



*£ ) w 2 ,Mio®# ) o4..A«S) 



Hence, the expression 1)240 represents the resumma- 
tion of the hi(Q 2 / $) logarithms over the intermediate 
Up scale, performed in such a way that both the loga- 
rithms In^fip/no) originating from the perturbative part 
of the DA and the ln(Q 2 / fi F ) logarithms from the hard- 
scattering part are resummed. The effect in the final 
prediction, at every order, is the same as if we had per- 
formed the complete renormalization-group resummation 
of ln(Q 2 //Xg) logarithms. 

Although by using the explicit results for J2J) and the 
evolution equation solution for <py, it is straightforward 
to employ 1)22(1 and obtain Th with resummed ln(fi F /Q 2 ) 
logs, it is much easier that the complete resummation 
is performed in the distribution amplitude. Hence, by 
adopting the common choice up = Q 2 , we avoid the 
need for the resummation of the ln(Q 2 /fj, F ) logarithms in 
the hard-scattering part, making the calculation simpler 
and hence, for practical purposes, the preferable form of 
F 17T (Q 2 ) is given by 



F in (Q 2 )=T H (x,Q 2 ,Q 2 ) ® $*(:r,Q 2 



(25) 



We stress here that in this approach, in which the con- 
sistent treatment of $ and Tjj dependence on (x F is re- 
quired, any other choice of /i 2 ? would lead to the same 
result, only the calculation would be more involved. 



V. CONCLUDING REMARKS 

We have sketched the higher-order PQCD calculational 
procedure for the hard exclusive quantities on the exam- 
ple of the pion transition form factor F 1 - K . 

Furthermore, we have argued that the Ky n prediction 
Q is independent of the factorization scale up at every 
order in as, when both the hard scattering part Th and 
the distribution amplitude $ are consistently treated re- 
garding the (i F dependence, i.e., in both quantities the 
In n F logarithms are resummed or in both quantities they 
are not resummed. The fip dependence of $ then exactly 
cancels the ptp dependence oiTn, the choice of the factor- 
ization scale is therefore nonessential and the predictions 
obtained by using any choice of fip are equal to the re- 
sults obtained using, for practical purposes the simplest 
intermediate choice up = Q 2 , where Q 2 represents the 
characteristic scale of the process. 

The true expansion parameter left is ag(/i^), with /j, 2 ^ 
representing the renormalization scale of the complete 
perturbatively calculable part of the pion transition form 
factor (|23|l . i.e., of 

T H {s,Q 2 ^ 2 )=T H {y,Q 2 ,Q 2 )®cpv{y,s,Q 2 , t i 2 ). (26) 

Therefore, although F 77! (Q 2 ) depends exclusively on the 
characteristic scale of the process Q 2 , we are left with 
the residual dependence on the /i 2 ^ scale, when calculat- 
ing to finite order. The intermediate scale at which the 
short- and long-distance dynamics separate, the factor- 
ization scale disappears from the final prediction at every 



order in as and therefore does not introduce any theoret- 
ical uncertainty into the PQCD calculation for exclusive 
processes. 

Above discussed calculational prescription for the fac- 
torization scale independent calculation is also upheld for 
other PQCD exclusive one-scale processes. However, in 
the case of exclusive processes which involve more than 
one typical scale the treatment of the factorization scale 
dependence is more involved. The subtlety in the pre- 
ceding approach lies in the fact that we have traded one 
dependence on fip for another one. Namely, the choice 
we have made there is that we resum In terms up to the 
relevant scale of the process Q 2 . Although, this might 
seem reasonable for the one external scale processes, such 
as the one which define F J7r (77* — > w) or the pion elec- 
tromagnetic form factor, in processes with two scales, for 
example in the one in which the general pion transition 
form factor Ry*^ (7*7* — ► ~k) appears, one immediately 



encounters the ambiguity of how to choose the relevant 
scale up to which the logs will be resummed. The exis- 
tence of such ambiguities seems to be unavoidable arte- 
fact of the PQCD calculation. 
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